Non-linear rheology of a nanoconfined simple fluid 
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We probe the rheology of the model liquid octamethylcyclotetrasiloxane (OMCTS) confined into 
molecularly thin films, using a unique Surface Forces Apparatus allowing to explore a large range 
of shear rates and confinement. We thus show that OMCTS under increasing confinement exhibits 
the viscosity enhancement and the non-linear flow properties characteristic of a sheared supercooled 
liquid approaching its glass transition. Besides, we study the drainage of confined OMCTS via the 
propagation of "squeeze-out" fronts. The hydrodynamic model proposed by Becker and Mugele 
[Phys. Rev. Lett. 91, 166104 (2003)] to describe such front dynamics leads to a conclusion in 
apparent contradiction with the dynamical slowdown evidenced by rheology measurements, which 
suggests that front propagation is not controlled by large scale flow in the confined films. 

PACS numbers: 64.70. pm, 83.50.Lh, 83.50-v, 83.60.Rs 



Supercooled liquids share qualitative rheological fea- 
tures upon approaching the glass transition [jy]: (i) their 
viscosity increases dramatically, and (ii) they exhibit 
non-newtonian properties (shear thinning) when the time 
scale of mechanical forcing becomes shorter than that of 
structural relaxation. The precise origin of such a be- 
havior is the subject of active theoretical and numerical 
investigations 0. Recently, an extension to How situa- 
tions of the mode-coupling theory (MCT) has been pro- 
posed, in order to describe this non-linear rheology Q. 
Now, a stringent test of theoretical predictions against 
experimental results requires measurements, over a large 
range of shear rate (7), of the non- linear properties as 
jamming is gradually approached. These are extremely 
challenging to perform on atomic glass formers, because 
of their elevated glass transition temperature and flow 
stress. The most comprehensive studies to date have fo- 
cused on colloidal suspensions of thermosensitive parti- 
cles, in which the volume fraction, hence the distance 
to jamming, can be finely tuned [J|. It has thus been 
shown that, in very good agreement with MCT, the flow 
stress of such suspensions exhibits a rate dependence all 
the weaker that the distance to glass transition is small, 
until a yield stress develops when the suspension gets 
jammed Such a behavior can be considered as the 
rheological hallmark of the approach to glass transition. 

Here, we show that increasing confinement represents 
an alternative pathway to bring a system close to its 
jammed state under well-controlled conditions. Surface 
Forces Apparatus (SPA) experiments have shown that 
simple liquids confined between solid walls below thick- 
nesses of a few molecular diameters exhibit enhanced 
flow resistance [H]. From SFA experiments probing the 
linear response of ultrathin liquid films, Demirel and 
Granick (DG) concluded to a confinement-induced dy- 
namical slowdown, akin to what occurs in supercooled 
liquids 6]. However, this conclusion has been challenged 
by other groups probing the large strain shear response 
of confined fluids [3, @|. Moreover, experiments by Becker 



and Mugele (BM) have shown that a confined liquid 
drains stepwise by expelling monolayers via the propa- 
gation of "squeeze-out fronts" Q. A model of the front 
dynamics, extending the work of Persson and Tosatti [lo| , 
led them to conclude that the confined fiuid retained its 
bulk viscosity, dissipation enhancement arising from high 
friction on the confining walls. 

The nature of the mechanisms by which the properties 
of liquids are affected by confinement at the molecular 
scale therefore remains an open question. Such an issue, 
which is of interest for the fundamental understanding of 
the jamming transition is also of paramount impor- 
tance for boundary lubrication and for nanofluidics, 
where the knowledge of the fiow properties of liquids con- 
fined into nanometer-sized channels or structures is cru- 
cial [13. 

In this Letter, we report on the first SFA study in 
which both large strain shear rheology and squeeze-out 
fronts measurements are performed, in the same exper- 
imental run, on the nonpolar liquid octamethylcyclote- 
trasiloxane (OMCTS), which has been used in the afore- 
mentioned works. 

(i) We show unambiguously, from fiow curves mea- 
sured over 6 decades of 7, that OMCTS under increasing 
confinement exhibits the viscosity enhancement and non- 
newtonian features of a supercooled liquid approaching 
the glass transition. 

(ii) We observe squeeze-out front dynamics in quanti- 
tative agreement with that previously reported [91 ■ When 
analyzed within the framework of the BM model, it re- 
sults in an effective viscosity two orders of magnitude 
lower than that directly measured in shear. We conclude 
that such an apparent contradiction arises from an im- 
proper assumption by BM about the nature of the mass 
transport mechanism at play during front propagation. 



Pxperiments were performed on a home-built SFA 14 1 
(Fig. [1]). The liquid is confined between two atomically 
smooth backsilvered mica sheets glued onto crossed cylin- 
drical lenses (radius of curvature ~1 cm). The normal 
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force Fn is measured by nieans of a load cell of stiffness 
9500 N.m"^. The "contact" area A, over which the mica 
sheets elastically flatten to form a circular parallel gap 
in which the liquid is confined, is monitored by videomi- 
croscopy. The thickness of the film, d, is determined by 
multiple beam interferometry [l^ and fast spectral cor- 
relation 14 , [l^ . Once confined under a given load, over 
an area A, the liquid is sheared, in a plane Couette ge- 
ometry, by moving laterally one surface at a velocity V 
in the range 10"'' — 10^ /im.s""', while measuring the re- 
sulting tangential force Ft with a cell of stiffness 5200 
N.m^^. The shear stress sensitivity of the instrument is 
Ft I A ~ 200 Pa. Our SFA has the unique feature of using 
the normal force signal as the input of a feedback loop, 
which allows to perform steady-state experiments over 
large shear amplitudes (up to hundreds of microns) un- 
der constant normal load conditions, whatever the level of 
confinement of the liquid. The mica sheets were prepared 
as described in , glued onto the cylindrical lenses us- 
ing a UV curing glue (NOA 81, Norland), and cleaved 
with adhesive tape immediately before being installed in 
the SFA, so as to obtain contaminant-free surfaces [isj . 
OMCTS from Fluka (purum grade >99%) was vacuum 
distilled before use. A drop (~ 150 /zL) of the hquid, fil- 
tered through a 0.2 /im membrane, was injected between 
the surfaces. It was then left at T = 20 ± 0.01°C for 12h 
in the sealed SFA, containing P2O5 to scavenge residual 
moisture, before beginning experiments. 
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FIG. 1: Force vs distance curve during approach of the sur- 
faces (loading velocity 0.5 nm.s~^). Inset: scheme of the 
setup. White light is shone on the confined film, and the 
transmitted intensity is sent (i) to a spectrometer for spectral 
analysis [l^, and (ii) to a CCD camera acquiring images of 
the contact area ^ at a rate of 55 s~^. 

Fig. [1] shows a force-distance profile measured upon 
quasi-static approach of the surfaces: it is clearly seen 
that below 6 nm[l^, the thickness of the confined liquid 
decreases by steps of approximately 8 A, which corre- 
sponds to the minor diameter of the slightly oblate OM- 
CTS molecule. This reflects the well-documented wall- 



induced layered structure of the fluid, which gives rise to 
the so-called solvation forces po| . 

We first focus on shear experiments performed on lay- 
ered OMCTS films with thicknesses ranging from 6 down 
to 2 monolayers. Over the whole range of confinement 
and velocity explored, we have observed: (i) a smooth 
stable shear response (see time trace in the inset of Fig. 
[21d), and (ii) a steady-state value of Ft which increases 
with V . On Fig. [2^, we plot the steady-state flow stress 
a = Ft I A versus shear rate — V jd for the different film 
thicknesses. The same data are plotted on Fig. [^b as the 
effective viscosity r^eff — <^/i versus 7. 
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FIG. 2: (a) a{j) for OMCTS films of (a) 6, (□) 5, (■) 4, (o) 3, 
and (•) 2 monolayers, (b) r;off(7), symbols as in (a). Insert: 
time trace of (•, left scale) a measured at V = O.l/im.s"^ 
on a 2nm-thick film, (line, right scale) the forth and back 
shear motion applied, (c) Master curve showing data from 
(b) plotted as rics/vo 7/7c. The solid hne is a fit of the 
form r)cs/rjd — 1/(1 + 7/7c)°'**. Inset: values for 770 (•, left 
scale) and l/7c (o, right scale) used for each n. 

It can be seen that, as the OMCTS thickness is reduced 
from n =6 to 2 monolayers: 
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(i) The flow stress, and hence the viscosity, steadily 
increases j2l| . 

(ii) The dependence of a on 7 shifts from hnear (new- 
tonian) to subhnear (shear-thinning). 

(iii) The crossover shear rate, 7c, above which non- 
newtonian behavior is observed, shifts to smaller values. 

(iv) For n <4 monolayers, power-law shear thinning 
(ct ~ 7", a < 1) at low 7 crosses over to a quasi-plateau 
regime. 

This set of features is characteristic of the approach to 
jamming, as reported experimentally |4| and predicted 
by MCT and numerical simulations [2, ISf . This is fur- 
ther supported by the fact that rjes('^) curves measured 
for different n collapse onto a single master curve when 
plotted as rjes/vo vs ■j/'jc, with 770 the zero shear viscos- 
ity (Fig. The inset of Fig. ^ shows that both 770 
and l/7c (i-e. the relaxation time of the liquid) sharply 
increase as the film thickness is decreased. The reduced 
viscosity obeys 7?eff/''7o — l/(l+7/7c)°'**, which is consis- 
tent with theoretical predictions for sheared supercooled 
systems Finally, the observation of a quasi-plateau 
regime which does not extend down to the lowest shear 
rates indicates that, in the present experiments, confined 
OMCTS approaches but does not reach jamming. This is 
consistent with the fact that (see inset of Fig. [5]d), upon 
cessation of shear, the stress relaxes (i) very slowly, over 
~ 5 s, and (ii) down to a non-measurable level. 

These observations lead us to conclude, in good agree- 
ment with DG 0, that OMCTS undergoes dynamical 
slowdown upon increasing confinement, similarly to a su- 
percooled system close above jamming. Such a conclu- 
sion contrasts with that of Klein Q or Israelachvili Q, 
who observed responses exhibiting a stick-slip dynam- 
ics which they interpret in terms of shear-melting of a 
confinement-induced ordered solid-like structure. Such a 
discrepancy might have two origins, (i) The use of differ- 
ent protocols for mica surface preparation: Indeed, the 
method employed in 0, Q , in contrast to that described 
above, may lead to surface contamination by a submono- 
layer of nanoparticles, which have been suggested as a 
possible reason for the observed stick-slip behavior [j^ . 
(ii) Differences in the crystallographic alignment of the 
confining surfaces, which is expected to affect the shear 
response of the intercalated molecular film [2^ . No sys- 
tematic investigation have been made so far of the effect 
of alignment between contaminant-free surfaces, and we 
therefore cannot discriminate between point (i) and (ii) 
above to explain differences. 

We now present the results from squeeze-out exper- 
iments. During loading, we record the light intensity 
transmitted through the contact area, along with F„ and 
d. We observe, as in [^, [l3|, that a film of thickness 
n monolayers drains via nucleation/growth of a circular 
region of thickness (n — 1) layers (see Fig. [3]). Nuclc- 
ation is accompanied by elastic relaxation of the confin- 
ing sheets, which are locally bent in the boundary zone 



connecting the regions of thickness n — 1 and n (Fig. [3] 
inset). This creates a 2D pressure gradient which then 
drives the monolayer expulsion The local curvature 
of the mica sheets induces a contrast in the transmit- 
ted intensity (Fig. [3^-c) which allows us to follow with 
time the position of the "squeeze-out" front. We have 
thus measured, for successive n — > n — 1 transitions, the 
squeeze-out time r needed to expel one monolayer from 
the contact area A. We have done so before and after 
rheology experiments, and did not observe any influence 
of shear history on front dynamics. 



10' 



(a) 


(b) 




(c) 



10 



14 



10' 



10 



12 



I I I I I I I I 



center 



•edgez 



boundary 



I I I I I I T I 



10^ 

10' 
10" 

10-' 



^ 10"- 
10"- 



4 
n 



Oh 











e) : 










c 






• 




5 





4 
n 



FIG. 3: (a-c): Sequence of images (96x96 /^m^) showing the 
front propagation during a 3 — > 2 transition, (d); /icB vs 
n (number of monolayers). (•) our results, (o) BM results, 
adapted from [gt]. Inset: schematic cross-section of the film 
during squeeze-out. (e): ?7eff vs n. (o) measured in shear, and 
(•) deduced from squeeze-out experiments. The horizontal 
line indicates the bulk viscosity of OMCTS. 

In the Persson and Tosatti (PT) model [10], the 
front velocity is related to the 2D pressure gradient by: 
Vp2D = -p2D/ieffT^, where p2D ^ Pa and p2D = pa 
[S (P — Fn/A is the applied pressure, p the fluid 
density and a the molecular size), V is the front velocity 
and ficS a viscous drag coefficient. The latter is deduced 
from the squeeze-out time as PcS ~ At:tP/ [pA). On 
Fig. [3Jl we have plotted /ies as a function of film thickness 
for our experiments, along with the values obtained by 
Becker and Mugele (BM) [9| . There is quantitative agree- 
ment between both data sets. BM have extended the PT 
model by assuming that front propagation is controlled 
by a "layered" Poiseuille flow between the front and the 
edge of the confinement area (Fig. [3]inset), and thus pro- 
posed that /Xcff should identify with the drag coefficient 
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of a Hele-Shaw flow, i.e. ^eS = 12^eff/(prf^), with fjos a 
shear viscosity and d the film thickness @. We use this 
expression to infer fj^ff(d) from the front dynamics. On 
Fig. [3^, we compare [2J] it to f7off((i) obtained from shear 
data. It appears that fjcs, which stays close to the bulk 
value down to 3-layer-thick films, is about two orders of 
magnitude lower than rj^ff. 

We propose the following explanation to this apparent 
paradox. Shear experiments are a straightforward way 
to measure rjes, in contrast to squeeze-out experiments, 
which require modelling of the front dynamics to infer 
a viscosity. Therefore, we consider that the reliable re- 
sults regarding rics{d) are those from shear rheology. We 
are then left with the observation of squeeze-out fronts 
which, given the ?7off(rf) obtained in shear, travel much 
faster than expected from the drag mechanism assumed 
in the BM picture, from which we have also drawn erro- 
neous conclusions in a recent study [17[. This suggests 
that the front dynamics is not controlled by the coher- 
ent sliding of adjacent incompressible molecular layers 
ahead of the front. Indeed, another piece of information 
emerges from the force-distance profile of Fig. [TJ between 
two steps, the film thickness is observed to decrease by 
about 3 A as the force is increased. Such a thickness 
variation is reversible upon load reduction. This shows 
that layered OMCTS films are substantially compress- 
ible, hence contain a non-negligible amount of free vol- 
ume, which is consistent with the fact that confined films 
do not reach jamming. This certainly facilitates local re- 
arrangements, and it is therefore likely that during prop- 
agation of a squeeze-out front, molecules in the layered 
region ahead of it permeate between layers in order to 
accommodate for density variations in the vicinity of the 
front. The apparent low resistance to front propagation 
suggests that permeation, rather than large scale coher- 
ent sliding of layers, controls mass transport ahead of the 
fronts. It implies that, pending further modelling, front 
dynamics cannot be used to infer a viscosity. 

In summary, we have probed the rheology of a simple 
fiuid under molecular confinement, and conclude that its 
behavior is akin to that of a sheared supercooled liquid 
close above the glass transition. This shows, as suggested 
by recent experiments on colloids that confinement 



which non-bulk behavior appears, and how it compares 
to the range of surface forces. 

We thank A. N. Morozov for fruitful discussion and C. 
Caroli for critical reading of the manuscript. 



can be used as an alternative route to finely control the 
approach to jamming. Our results now raise two im- 
portant questions, (i) We observe a liquidlike behavior 
down to the thinnest film investigated, which brings up 
the issue of how to cross the jamming transition under 
confinement. Two routes can be envisaged. It can be 
done by varying the chemical corrugation of the walls, as 
shown in friction experiments 25| or in numerical sim- 
ulations [26], or, as mentioned above, by changing the 
orientation between the crystalline lattices of the confin- 
ing surfaces, (ii) We find that 6 layer-thick films already 
exhibit a viscosity two orders of magnitude larger than 
the bulk value. This raises the question of the scale below 
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